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A search for nonresonant new phenomena, originating from either contact interactions or large 
extra spatial dimensions, has been carried out using events with two isolated electrons or muons. 
These events, produced at the LHC in proton-proton collisions at y/s = 7 TeV, were recorded by the 
ATLAS detector. The data sample, collected throughout 201 1 , corresponds to an integrated luminosity 
of 4.9 and 5.0 fb _1 in the e+e~ and channels, respectively. No significant deviations from the 
Standard Model expectation are observed. Using a Bayesian approach, 95% confidence level lower 
limits ranging from 9.0 to 13.9 TeV are placed on the energy scale of Uqq contact interactions in the 
left-left isoscalar model. Lower limits ranging from 2.4 to 3.9 TeV are also set on the string scale in 
large extra dimension models. After combining these limits with results from a similar search in the 
diphoton channel, slightly more stringent limits are obtained. 
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A search for nonresonant new phenomena, originating from either contact interactions or large 
extra spatial dimensions, has been carried out using events with two isolated electrons or muons. 
These events, produced at the LHC in proton-proton collisions at *Js = 7 TeV, were recorded by 
the ATLAS detector. The data sample, collected throughout 2011, corresponds to an integrated 
luminosity of 4.9 and 5.0 fb _1 in the e + e~ and fi + fi~ channels, respectively. No significant deviations 
from the Standard Model expectation are observed. Using a Bayesian approach, 95% confidence level 
lower limits ranging from 9.0 to 13.9 TeV are placed on the energy scale of Uqq contact interactions 
in the left-left isoscalar model. Lower limits ranging from 2.4 to 3.9 TeV are also set on the string 
scale in large extra dimension models. After combination of these limits with results from a similar 
search in the diphoton channel, slightly more stringent limits are obtained. 
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I. INTRODUCTION 

Extensions to the Standard Model (SM), such as 
quark/lepton compositeness and large extra dimensions, 
predict modifications to the SM dilepton invariant mass 
spectra. This paper presents a comparison of the num- 
ber of expected and observed events at high mass in the 
dielectron and dimuon datasets collected by the ATLAS 
detector [l[ in 2011. These events resulted from proton- 
proton collisions produced at i/s = 7 TeV by the LHC Q ■ 
The data are interpreted in the context of contact inter- 
actions (CI) and virtual graviton exchange in the Arkani- 
Hamed-Dimopoulos-Dvali (ADD) model Q. 

In the SM, quarks and leptons are fundamental par- 
ticles. However, if they are composite particles, with at 
least one common constituent, the interactions of these 
constituents would manifest themselves through an ef- 
fective four-fermion contact interaction at energies well 
below the compositeness scale. This type of contact inter- 
action could also describe a new force with a messenger 
too heavy for direct observation at the LHC, in analogy 
with Fermi's nuclear fj decay theory |4(. 

The Lagrangian for a general contact interaction has 
the form [5j: 

+ »7rr ^rTm^R^rT'Vr C 1 ) 
+2j7lr VlTmV'l tp-BTf^ipR } , 

where g is a coupling constant chosen so that g 2 /47r = 1; 
A is the contact interaction scale, which in the context 
of compositeness models, is the energy scale below which 
fermion constituents are bound; and ipL.R are left-handed 
and right-handed fermion fields, respectively. The pa- 
rameters rjij, where i and j are L or R (left or right), 
define the chiral structure of the new interaction. Spe- 
cific models are constructed by assigning particular com- 
binations of these parameters to be — 1, or +1. For 
example, the left-left isoscalar model (LLIM) is defined 
by setting ?7 LL = ±1 and ?7rr = ?7lr = 0. The LLIM 



model, commonly used as a benchmark for contact inter- 
action searches [6J, is utilized in this analysis. 

The addition of the contact interaction Lagrangian to 
that of the SM modifies the Drell-Yan (DY) production 
cross section (qq — > Z/j* — > £ + £~). The largest devia- 
tions in the dilepton invariant mass spectra, either con- 
structive or destructive, are expected at high mass and 
are determined by the sign of the parameter rjij and the 
scale A. The differential cross section for the process 
qq — > £ + £~ , including a contact interaction, can be sep- 
arated into three components: a SM DY term, a pure 
contact interaction term (Fq) and a DY-CI interference 
(Fi) term: 

dcr _ dopy _ Fi{m u ) F c (m u ) 
dm u ~ dm u ^ A 2 A 4 ' U 

where m« represents the final-state dilepton mass. The 
full form of this expression is given in Ref . 0] ■ Construc- 
tive (destructive) interference corresponds to ?7ll = — 1 
(+1). At the largest A values to which this analysis is 
sensitive, both interference and pure contact interaction 
terms play significant roles. For example, at dilepton 
masses greater than 400 GeV and A = 12 TeV, the mag- 
nitude of the middle term in Eq. ^ , which depends on 
the interference, is about twice that of the last term. 

Nonresonant deviations in the high mass dilepton in- 
variant mass spectra are also predicted in large extra 
dimension models. These models were introduced to 
address some of the major unresolved issues in parti- 
cle physics such as the hierarchy problem. The latter 
deals with the question of why gravity appears weak in 
comparison to the other three SM interactions and why 
the electroweak scale (~1 TeV) is 16 orders of magni- 
tude smaller than the Planck scale (Afpi ~ 10 16 TeV). 
Arkani-Hamed, Dimopoulos, and Dvali addressed these 
issues by postulating the existence of n flat additional 
spatial dimensions of common size R, compactified on an 
n-dimensional torus The fundamental Planck scale 
in (4+n)-dimensional spacetime, Md, is then related to 
the scale M Pl by Gauss's law: Af|, = M£ +2 R n . Con- 
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sequently, the hierarchy problem can be solved with a 
Planck scale resulting from a fundamental scale (Mg ) 
near 1 TeV if the volume (oc R n ) is large enough. 

In the ADD model, the SM particles and their inter- 
actions are confined to a three-dimensional slice of the 
multidimensional world, but gravity permeates the addi- 
tional dimensions of size R. This results in Kaluza-Klein 
(KK) modes of the graviton. The mass splitting of these 
KK modes is determined by the factor 1/R. Resolution 
of the hierarchy problem necessitates large extra dimen- 
sional volumes and consequently implies small values of 
1/R. This results in an almost continuous spectrum of 
KK graviton states and hence a nonresonant increase in 
the expected rate of dilepton events at large invariant 
mass. Performing the sum over the KK modes in the vir- 
tual graviton exchange process leads to an integral which 
has to be regulated by an ultraviolet cutoff value (Atjv)- 
The ADD model is a low-energy effective theory valid 
below the scale of the onset of quantum gravity, charac- 
terized by the scale Ms . The convention used throughout 
this analysis is to equate the cutoff to the scale of the ef- 
fective theory (Atjv=Ms). 

For virtual graviton exchange, it is standard practice 
to present limits on the size of the extra dimensions in 
terms of Ms , taken to be the string scale, which is related 
to Md by the following expression |8j: 



Ms 



2y^ 



1/0+2) 



M D 



(3) 



The strength of gravity in the presence of extra dimen- 
sions is typically parametrized by ?7g = J 7 / Mg, where T 
is a dimensionless parameter of order unity. The defini- 
tion of T depends on the formalism chosen [8[ , with three 
popular conventions: Giudice-Rattazzi- Wells (GRW) 
Hewett [13 and Han-Lykken- Zhang (HLZ) [nj]. The dif- 
ferent values are 



T 
T 

T 



1, (GRW) 
2A ±2 



71-2 



(Hewett) 

71 

for n > 2. (HLZ) 



(4) 



In the GRW and HLZ representations, gravitational 
effects interfere constructively with the SM processes, 
while in Hewett 's convention there can be destructive 
or constructive interference. This is encapsulated in the 
parameter A, which is equal to +1 (—1) for constructive 
(destructive) interference. 

The total cross section (cr to t), including effects of 
qq- and gg-initiated virtual graviton exchange, may be 
parametrized as 



Ctot = °SM + m F int +VgFg, 



(5) 



where ctsm is the SM cross section for the process being 
considered, and i<l n t and Fq are functions of the cross sec- 
tions involving the interference and pure graviton effects, 
respectively. Note that the interference term has a linear 



dependence on t/q, and therefore a quartic dependence 
on Ms (i.e., ?yc oc 1/Mg), whereas the pure graviton ex- 
change term is quadratic in r/Q and therefore has a 1/M| 
dependence. A study of signal yields in the kinematic 
range relevant to this analysis shows that the pure gravi- 
ton term dominates those yields. This is in part due to 
the fact that the gg-initiated contribution to the graviton 
exchange process does not interfere with the gq-initiated 
DY process. Results are nevertheless presented for both 
1/M| and 1/Mf priors. 

Previous searches for contact interactions have been 
carried out in neutrino-nucleus and electron-electron 
scattering [l2|, EH; as wen as at electron-positron |l4- 
181. electron- proton [l9|, [2(| , and hadron colliders [21- 
28j . In the case of eeqq contact interactions, the limits 



in the LLIM for all quark flavors from e + e experiments 
are A" > 7.2 TeV and A+ > 12.9 TeV [3 at 95% con- 
fidence level (C.L.) for 7/ll = —1 and +1, respectively. 
These limits assume that contact interactions of electrons 
with all quark flavors are of the same strength. The best 
limits set in the specific case of first generation quarks 
are A~ > 9.1 TeV and A+ > 8.6 TeV QJ] at 95% C.L. 
In the case of eeqq contact interactions, the best limit 
for constructive interference is A~ > 10.1 TeV from the 
ATLAS analysis of the first 1 fb" 1 of 2011 data H. The 
best limits in the case of fifiqq contact interactions are 
from an analysis of the same data: A - > 8.0 TeV and 
A+ > 7.0 TeV [13. 

Previous searches for large extra dimensions in the 
ADD model via virtual graviton exchange have been per- 
formed at electron-positron I29ll33 . electron-proton [2(| 
H3 , and hadron colliders [25L l36M42j | . Presently, the most 
stringent mass limits in the dielectron and dimuon chan- 
nels require Ms > 2.8 TeV for each channel and 3.1 TeV 
when combined (in the GRW formalism with no K fac- 
tor) [§3- The best limits to date arise from the combi- 
nation of these dilepton results with those from a search 
in the diphoton final state, which increases the limit by 
~ 0.1 TeV [38]. The following sections describe the first 
virtual graviton exchange search performed by ATLAS 
using dilepton data and its combination with an ATLAS 
diphoton data search [42| . 



II. ATLAS DETECTOR 

ATLAS is a multipurpose particle detector composed 
of three main subsystems: the inner tracking detector, 
the calorimeter system and the muon spectrometer. The 
inner detector is used to track charged particles within a 
pseudorapidity ?£] in the range \rj\ < 2.5. It comprises a 



1 ATLAS uses a right-handed coordinate system with its origin at 
the nominal interaction point (IP) in the center of the detector 
and the 2-axis along the beam pipe. The rr-axis points from the 
IP to the center of the LHC ring, and the j/-axis points upward. 
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silicon pixel detector, a silicon-strip tracker and a tran- 
sition radiation tracker. An axial 2 T magnetic field is 
produced by a superconducting solenoid. 

The calorimeter system, covering \rj\ < 4.9, surrounds 
the solenoid and provides three-dimensional reconstruc- 
tion of electromagnetic and hadronic showers. The 
lead/liquid-argon electromagnetic sampling calorimeter 
covers \r/\ < 2.5 and is finely segmented with a readout 
granularity varying by layer and with cells as small as 
0.025 x 0.025 in (ij, <j)) to provide precise energy and po- 
sition resolution, as needed for electron and photon iden- 
tification and energy measurement. Hadron calorimctry 
is provided by an iron/scintillator tile calorimeter in the 
central pseudorapidity range \rj\ < 1.7 and a lead/liquid- 
argon calorimeter extending the pseudorapidity range up 
to |f?| = 3.2. Both the electromagnetic and hadronic 
calorimeters have liquid-argon-based forward detectors, 
with copper or tungsten as an absorber, to extend cover- 
age up to |?7| = 4.9. 

Outermost is the muon spectrometer, another key de- 
tector component for this analysis. Three layers of preci- 
sion tracking chambers, comprising monitored drift tubes 
and cathode strip chambers, enable muon reconstruction 
up to = 2.7. The magnetic field is provided by three 
large air-core toroidal magnet systems (one barrel and 
two end caps), each consisting of eight azimuthally sym- 
metric superconducting coils. Triggering capability up to 
= 2.4 is provided by fast resistive plate chambers in 
the barrel and thin-gap chambers in the end caps. 



III. SIGNAL AND BACKGROUND MODELING 

The dominant background contribution comes from 
the SM DY process with smaller contributions from tt 
and electroweak diboson (WW, WZ, and ZZ) produc- 
tion. In the dielectron channel, there is also a signif- 
icant background from multijet and W+jets events in 
which jets are misidentified as electrons. Backgrounds 
are estimated using fully simulated Monte Carlo (MC) 
samples except for the combined multijet and W+jets 
background, which is determined from the data. 

DY samples are generated with Pythia 6.421 [43[ 
using MRST2007 LO** parton distribution functions 
(PDFs) 4J|. The diboson background is generated with 
Herwig 6.510 [H using MRST2007 LO** PDFs. For 
the tt background, event generation is performed with 
Mc@nlo 4.01 [46] and the CTEQ 6.6 PDFs [47|, as well 
as Herwig to model the underlying event and parton 
showers. Production of diboson and tt events relies on 
Jimmy 4.31 [48[ to describe multiple interactions. 

For the contact interactions analysis, Pythia 6.421 
and the MRST2007 LO** PDFs are used to generate 



Cylindrical coordinates (r, <j>) are used in the transverse plane, <j> 
being the azimuthal angle around the beam pipe. The pseudora- 
pidity is defined in terms of the polar angle 8 as r\ = — lntan(0/2). 



both DY and CI as part of the same process in order 
to properly handle the interference between DY and CI 
contributions, as shown in Eq. @. All quark flavors 
contribute to the DY+CI production. 

Final-state radiation is simulated with Photos [4£| 
for all the samples listed above. Higher-order QCD cor- 
rections are implemented via mass-dependent K factors 
defined as the ratio of the next-to-next-to-leading-order 
(NNLO) Z/7* cross section calculated using Phozpr 
and the MSTW2008 NNLO PDFs [H|, to the LO Z '/ 7 * 
cross section from Pythia. Higher-order electroweak 
corrections originating from virtual gauge-boson loops 
are computed using the Horace NLO event genera- 
tor (52j. The mass-dependent QCD and electroweak K 
factors are applied to both DY and DY+CI samples. 

For the large extra dimensions analysis, Sherpa 
1.3.1 [H[ and the CTEQ6L PDFs H are used to gen- 
erate DY+ADD events at leading order. The generated 
dilepton mass is required to be less than the scale Ms, 
which is chosen to be in the range between 2 and 3 TeV in 
this study, since the model is not valid at energies beyond 
that scale. 

Finally, the generated samples are processed through 
a full simulation of the ATLAS detector 55] based on 
Geant 4 [56[ and reconstructed with the same software 
as for the data. Several corrections derived from data 
control samples are applied to the simulated samples. 
Specifically, the energy scale and resolution for electrons 
are corrected so that the Z — > ee mass distribution in 
simulation matches the data [57j ] . Similarly, the muon 
transverse momentum scale and resolution are adjusted 
to reproduce the muon tracking performance measured 
in Z — y \i\x data as well as several dedicated data sam- 
ples [HI . The effect of pileup (multiple pp interactions in 
the same or nearby bunch crossings) is included by super- 
imposing minimum bias events, in the same proportion 
as in data, on top of the hard scattering process gener- 
ated. Small corrections are included in the analysis to 
properly describe the pileup conditions for the selected 
data sample. 



IV. EVENT SELECTION AND COMPARISON 
BETWEEN EXPECTED AND OBSERVED 
YIELDS 

This analysis follows the same event selection as the 
search for new heavy resonances [5^] in the dielectron 
channel but uses a somewhat tighter selection in the 
dimuon channel. 

The data sample was collected during LHC operation 
in 2011 and corresponds to a total integrated luminosity 
of 4.9 and 5.0 fb _1 in the e + e~ and /j, 4 final states, re- 
spectively. The events recorded by the ATLAS detector 
were selected by requiring that they pass specific trig- 
gers. The trigger for the dielectron dataset required the 
presence of two electromagnetic clusters consistent with 
originating from electrons with transverse momentum 
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above 20 GeV, whereas events in the muon dataset were 
required to pass at least one of two single-muon triggers 
with p T thresholds of 22 GeV and 40 GeV. 

After passing the trigger selection, events are required 
to have a pair of either electrons or muons with px greater 
than 25 GeV. Furthermore, events are required to be 
recorded during stable beam conditions and with detec- 
tor components operational. To reject cosmic ray events 
and beam halo background, events are required to have 
a reconstructed vertex with at least three charged par- 
ticle tracks with pt > 0.4 GeV. If more than one such 
vertex is found, the vertex with the largest £p T is se- 
lected as the primary vertex of the event, where the sum 
is over all charged particles associated with the given ver- 
tex. Electron candidates are confined to \rj\ < 2.47, with 
the calorimeter barrel-to-end-cap transition region 1.37 
< \rj\ < 1.52 excluded due to the degraded energy resolu- 
tion in this region. No explicit rj requirement is placed on 
muon candidates, but the selection described below leads 
to negligible acceptance beyond \r/\ of approximately 2.5. 

Electron candidates are formed from clusters of cells 
in the electromagnetic calorimeter where energy is de- 
posited. Identification criteria based on the transverse 
shower shape, the leakage into the hadronic calorime- 
ter, and the association to an inner detector track are 
applied to the cluster to satisfy the medium electron def- 
inition 57]. The electron energy is obtained from the 
calorimeter measurements and its direction from the as- 
sociated inner detector track. A hit in the first layer of 
the pixel detector is required (if an active pixel module 
is traversed) to suppress background from photon con- 
versions. Further jet background suppression is achieved 
by demanding that the highest-px electron in the event 
be isolated. To this effect, the sum of the transverse 
energies, S-Eti in calorimeter cells within a radius R = 
y/ (Arj) 2 + (A<j)) 2 of 0.2 around the electron direction, is 
required to be less than 7 GeV. The core of the electron 
energy deposition is excluded and the sum is corrected 
for transverse shower leakage and pileup. The two elec- 
tron candidates are not required to have opposite charge 
because of possible charge misidentification either due to 
bremsstrahlung or limited momentum resolution of the 
inner detector at high pt- If the event contains more 
than two selected electrons, the two electrons with the 
highest-px sum are chosen. For these selection criteria, 
the overall event acceptance for DY events has a small 
dependence on the dielectron mass above 500 GeV, with 
a value of approximately 65% at 1 TeV. 

Muon candidates are reconstructed independently in 
the inner detector and the muon spectrometer. The mo- 
mentum is taken from a combined fit to the measure- 
ments from the two subsystems. To obtain optimal mo- 
mentum resolution and accurate modeling by the simula- 
tion, muon candidates are required to have at least three 
hits in each of the inner, middle, and outer detector layers 
of the muon spectrometer, and to have at least one hit in 
each of two different layers in the nonbending xy plane. 
To suppress background from cosmic rays, requirements 



are imposed on the primary vertex (PV) position and 
the muon impact parameter relative to the PV: z co- 
ordinate of the PV |zpy| < 200 mm, muon transverse 
impact parameter \do\ < 0.2 mm and muon z coordi- 
nate \z — zpv| < 1 mm. Furthermore, the muons are 
required to be isolated to reduce background from jets: 
Xpt(-R < 0-3)/pt(m) < 0.05, where the sum is over inner 
detector tracks within a radius of 0.3 around the muon 
direction. If more than one opposite-sign muon pair is 
found in an event, the pair with the highest-px sum is 
chosen. The overall event acceptance for DY events has 
only a weak dependence on the dimuon mass, with a 
value of approximately 40% at 1 TeV. This is lower than 
the acceptance in the dielectron channel primarily due to 
the stringent requirements on the presence of hits in all 
three layers of the muon spectrometer and the extent of 
the three-layer geometrical coverage. 

The Vy+jets background in the dimuon channel is esti- 
mated from simulated samples and is found to be negligi- 
ble since the event must contain two well-measured high- 
Pt isolated muons. Likewise, the multijet background, 
estimated directly from the data by reversing the muon 
isolation criterion, is found to be negligible. The mul- 
tijet and W+jets backgrounds are not negligible in the 
dielectron channel. They are estimated primarily from 
the data using several methods [591 ] ■ The first method 
determines the multijet background from the data and 
relies on the MC simulation for the VF+jets contribu- 
tion. The background is measured with a template built 
by reversing one of the electron identification criteria and 
normalized to data in the range 70 < m ee < 200 GeV. 
Another independent method that is sensitive to both 
multijet and VF+jets backgrounds uses jet-enriched data 
samples either from jet triggers or from the same trigger 
used to select the events in this analysis. The method 
relies on jet misidentification rates, defined as the num- 
ber of jets that pass the full electron selection divided by 
the number that pass a loose electron selection obtained 
by reversing one of the identification criteria. The back- 
ground estimate is then constrained by a fit in the range 
140 < m ee < 850 GeV. The final combined multijet and 
Vy+jets background is obtained with a simple average of 
the expected event yields from the different methods. 

Extensive comparisons between data and MC simula- 
tion were performed at the level of single-lepton distri- 
butions to confirm that the simulation reproduces the 
selected data, especially at high momentum. Figure Q] 
shows good data-MC agreement in the lepton transverse 
momentum distributions for events passing all selection 
criteria. 

Figure [2] shows the dielectron and dimuon mass distri- 
butions for selected events. Also shown are the predicted 
contributions from SM and new phenomena (NP) for sev- 
eral choices of model parameters. The expected SM dis- 
tribution is dominated by the DY process over the entire 
mass range and is found to describe the data well. The 
level of agreement with the SM expectation is also illus- 
trated in Fig. [31 which shows the number of events above 
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Electron p T [GeV] 




FIG. 1. Lepton transverse momentum distributions in the 
dielectron (top panel) and dimuon (bottom panel) channels 
for data (points) and Monte Carlo simulation (histograms). 
The bin width is constant in log(pT). 



a minimum mass m™ n . 

A more quantitative comparison is provided in Ta- 
bles U and |TI] showing the numbers of observed and ex- 
pected events in the dielectron and dimuon channels, re- 
spectively. The expected yields are normalized to the 
number of events observed in the Z peak control region 
(70 < mu < 110 GeV). The mass region shown in these 
tables corresponds to the CI search region defined by 
mu > 400 GeV. These tables also display the expected 
yields for the SM+CI signal for the two scenarios where 
the CI interferes either constructively or destructively 
with the SM. 



V. SYSTEMATIC UNCERTAINTIES 

Except for the multijet and VF+jets background con- 
tributions to the dielectron channel, all signal and back- 
ground event yield estimates are based on MC simulation. 
Because these yields are normalized in the Z peak control 
region, only mass-dependent systematic uncertainties af- 
fect the event yield estimates in the higher-mass signal 




80 100 200 300 400 1000 2000 3000 

m ee [GeV] 




80 100 200 300 400 1000 2000 3000 



rrv [GeV] 

FIG. 2. Dielectron (top panel) and dimuon (bottom panel) in- 
variant mass distributions for data (points) and Monte Carlo 
simulation (filled histograms). The open histograms corre- 
spond to the distributions expected in the presence of con- 
tact interactions or large extra dimensions for several model 
parameters. The bin width is constant in log(m«). 



region. The only exception is a 5% uncertainty applied 
to the signal yield to account for the uncertainty in the 
Z/7* cross section which affects the signal normalization. 

Experimental uncertainties arise from lepton en- 
ergy/momentum scale and resolution, as well as trigger, 
reconstruction and identification efficiencies. In the di- 
electron channel, the largest experimental uncertainty 
comes from the combined multijet and W^+jets back- 
ground estimate. It is determined from the envelope 
of the three separate methods used, including the effect 
of varying the mass ranges in the background fits and 
the uncertainties in the 77 and pt dependence of the jet 
misidentification rates. Electron energy scale and resolu- 
tion are determined from data via J/ip — > ee and Z — » ee 
mass distributions, as well as studies of electron E/p in 
W — >• ev decays [53]. The uncertainty in the constant 
term that dominates the resolution at high energy has 
negligible impact on this analysis. A somewhat larger 
impact comes from the energy scale knowledge, result- 
ing in a systematic error of 1.2% and 2.4% for dielectron 
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TABLE I. Expected and observed numbers of events in the dielectron channel for the contact interactions search region. The 
yields are normalized to the Z peak control region and include predictions for SM backgrounds as well as for SM+CI with 
different CI scales for constructive (A - ) and destructive (A + ) interference. The errors quoted originate from both systematic 
uncertainties and limited MC statistics. 



m ee [GeV] 


400-550 


550-800 


800-1200 


1200-1800 


>1800 


DY 


203 ± 10 


62.5 ± 3.4 


12.1 ± 0.9 


1.38 ± 0.17 


0.085 ± 0.025 


ti 


22.6 ± 2.1 


4.05 ± 0.34 0.308 ± 0.026 


< 0.05 


< 0.01 


Diboson 


12.1 ± 0.7 


4.08 ± 0.21 


0.88 ± 0.05 


0.111 ± 0.006 


0.0100 ± 0.0006 


Multijet/W+jets 


38 ± 23 


11 ± 8 


2.0 ± 1.8 


0.24 ± 0.28 


0.022 ± 0.029 


Total background 


276 ± 25 


82 ± 9 


15.3 ± 2.0 


1.74 ± 0.33 


0.12 ± 0.04 


A~ = 3 TeV 


1460 + 70 


1400 ± 80 


1090 ± 60 


525 + 35 

'J t-> 'J \ 'J' J 


148 ± 13 


A~ = 4 TeV 


680 + 40 


519 ± 27 


360 ± 21 


171 ± 12 


44 ± 4 


A~ = 5 TeV 


463 ± 30 


281 ± 15 


162 ± 9 


77 ± 5 


19.8 ± 1.9 


\- — 7 TeV 


339 + 97 


145 ± 10 


59 ± 4 


99 + 1 6 

_ _ . 1 / 1 -L.LI 


4.8 ± 0.5 


A^ = 12 TeV 


293 ± 27 


96 ± 9 


23.6 ± 2.3 


5.1 ± 0.5 


0.87 ± 0.14 


A+ = 3 TeV 


1080 ± 50 


1120 ± 60 


920 ± 50 


493 ± 33 


128 ± 11 


A+ = 4 TeV 


484 ± 30 


373 ± 20 


291 ± 17 


156 ± 10 


40 ± 4 


A+ = 5 TeV 


342 ± 27 


182 ± 11 


114 ± 6 


61 ± 4 


18.3 ± 1.6 


A+ = 7 TeV 


268 ± 27 


102 ± 10 


37.4 ± 2.6 


15.1 ± 1.0 


4.3 ± 0.4 


A+ = 12 TeV 


260 ± 27 


82 ± 9 


15.1 ± 2.2 


2.5 ± 0.4 


0.41 ± 0.08 


Data 


270 


88 


17 


3 






TABLE II. Expected and observed numbers of events in the dimuon channel for the contact interactions search region. The 
yields are normalized to the Z peak control region and include predictions for SM backgrounds as well as for SM+CI with 
different CI scales for constructive (A - ) and destructive (A + ) interference. The errors quoted originate from both systematic 
uncertainties and limited MC statistics. 



m MM [GeV] 


400-550 


550-800 


800-1200 


1200-1800 


>1800 


DY 


123 ± 6 


37.4 ± 2.2 


7.1 ± 0.6 


0.82 ± 0.11 


0.058 ± 0.022 


ti 


13.4 ± 1.4 


3.1 ± 0.5 


0.04 ± 0.12 


< 0.05 


< 0.01 


Diboson 


7.9 ± 0.4 2.66 ± 0.15 0.55 ± 0.04 0.075 ± 0.006 0.0124 ± 0.0031 


Total background 


145 ± 6 


43.2 ± 2.2 


7.7 ± 0.6 


0.89 ± 0.11 


0.070 ± 0.022 


A^ = 3 TeV 


870 ± 50 


770 ± 50 


580 ± 40 


296 ± 28 


82 ± 22 


= 4 TeV 


405 ± 19 


301 ± 17 


201 ± 14 


87 ± 8 


27 ± 7 


A" = 5 TeV 


256 ± 12 


159 ± 8 


94 ± 6 


41 ± 4 


12.7 ± 3.4 


A" = 7 TeV 


184 ± 9 


79 ± 4 


30.1 ± 1.9 


12.3 ± 1.2 


2.9 ± 0.8 


A" = 12 TeV 


157 ± 9 


50.6 ± 3.1 


12.3 ± 0.9 


2.81 ± 0.31 


0.53 ± 0.15 


A+ = 3 TeV 


628 ± 31 


650 ± 40 


500 ± 40 


248 ± 23 


75 ± 20 


A+ = 4 TeV 


271 ± 12 


203 ± 11 


159 ± 11 


85 ± 8 


22 ± 6 


A+ = 5 TeV 


182 ± 9 


98 ± 5 


64 ± 4 


31.4 ± 2.9 


11.5 ± 3.0 


A+ = 7 TeV 


141 ± 8 


50.8 ± 3.1 


19.7 ± 1.2 


8.4 ± 0.8 


2.5 ± 0.7 


A+ = 12 TeV 


140 ± 8 


40.2 ± 3.0 


7.4 ± 0.7 


1.57 ± 0.20 


0.25 ± 0.08 


Data 


151 


36 


9 


1 






masses of 1 and 2 TeV, respectively. A slight efficiency 
drop of 1.0% per TeV is predicted by the simulation due 
to the isolation requirement on the leading electron. To 
account for this effect, an uncertainty of the same mag- 
nitude is introduced. 



In the dimuon channel, the largest contribution to the 
experimental systematic error comes from the muon re- 
construction efficiency and muon resolution. A slight 
drop in reconstruction efficiency is predicted by the simu- 
lation at high pt due to the presence of additional hits in 
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FIG. 3. Distribution of the number of events with dilepton 
mass above mJJ 1 " for data (points) and SM prediction from 
Monte Carlo simulation (filled histograms, shaded gray) in the 
dielectron channel (top panel) and dimuon channel (bottom 
panel). The open solid and dashed histograms correspond to 
the expected distributions in the presence of contact interac- 
tions or large extra dimensions for several model parameters. 
The bin width is constant in log(mJ3 m ). 



the muon spectrometer from muons undergoing large en- 
ergy loss in the detector. An uncertainty of 3.0% (6.0%) 
at a dimuon mass of 1 (2) TeV is assessed, corresponding 
to the magnitude of this effect. The limited knowledge 
of the momentum scale determined from Z — > data 
has a negligible impact on the analysis. The momen- 
tum resolution in the simulation is adjusted based on 
Z — >• /i/i and W — > [iv data, as well as dedicated straight 
muon track data collected with the toroids turned off and 
tracks passing through overlapping sectors in the muon 
spectrometer. The latter provide two independent mo- 
mentum measurements for the same muon. The toroid- 
off and overlapping sector tracks are key to determining 
the muon reconstruction performance at high px- The 
uncertainty in the muon resolution, taken as equal in 
magnitude to the correction applied to the simulation, 
results in a change in the event yield of 1.2% (12%) for 
m m = 1 (2) TeV. 



The largest error contribution due to theory arises from 
limited knowledge of the PDFs, as, and QCD K factors. 
Scale uncertainties are computed by taking the maximum 
deviations obtained by independently varying the renor- 
malization (/ir) and factorization (/xf) scales by a factor 
of 2 but with the constraint that the ratio (Uf//Ur does 
not change by more than a factor of 2. The as and PDF 
uncertainties are determined with the MSTW2008NNLO 
eigenvector PDF sets and the different PDFs correspond- 
ing to variations of as. The overall uncertainty is com- 
puted using 90% C.L. ranges and includes the enve- 
lope of the uncertainty bands for the following different 
PDF sets: MSTW2008, NNPDF2.1, CT10, and CT10W. 
PDFs are the largest source of uncertainty, with the enve- 
lope of all PDFs considered becoming the dominant con- 
tribution above a few hundred GeV. Uncertainties in the 
electroweak K factor [6(| originate from the calculation 
of real boson radiation, O(aas) corrections, higher-order 
electroweak corrections, an assumed uncertainty of 10% 
in the contribution from photon-induced processes, and 
a difference in the electroweak renormalization scheme 
definition used in Pythia and in the calculation of the 
electroweak corrections with Horace. The latter source 
is the largest contribution to the electroweak uncertainty. 

The systematic uncertainties are summarized in Ta- 
ble IIIII Although not explicitly listed in this table, the 
uncertainty due to limited MC statistics is also taken 
into account in the limit setting. For DY+CI MC sam- 
ples, this uncertainty grows from about 4% at low m« to 
about 30% at high m u for A = 12 TeV. 



VI. STATISTICAL ANALYSIS 

The data analysis proceeds with a Bayesian method 
to compare the observed event yields with the expected 
yields for a range of different NP model parameters 
(where the NP corresponds to either contact interactions 
or large extra dimensions). Specifically, the number of 
expected events in a given search region is 

A* = ™DY+Np(#, V) + n non -DY bg(^) , (6) 

where ?i.dy+np($, v) is the number of events predicted 
by the DY+NP simulation for a particular choice of NP 
model parameter 9, n non _DY bg(^) is the number of non- 
DY background events, and v represents the set of Gaus- 
sian nuisance parameters that account for systematic un- 
certainties. The parameter 9 corresponds to a choice of 
energy scale A and interference parameter ?7ll in the CI 
analysis or to a choice of string scale Ms and formalism 
in the ADD analysis. In the case of the CI analysis, the 
input to evaluate the complete set of /i values is shown in 
Tables Q] and [TT] for the dielectron and dimuon channels, 
respectively. For each mass bin, a second-order polyno- 
mial is used to model the dependence of /i on 1/A 2 . In 
the case of the ADD analysis, [i is also parameterized by 
a second-order polynomial but as a function of 1/M|. 
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TABLE III. Summary of systematic uncertainties in the expected numbers of events for a dilepton mass of 1 TeV (2 TeV). NA 
indicates that the uncertainty is not applicable. 



Source 




ee 




/'/' 




Signal 


Back 


ground 


Signal 


Background 


Normalization 


5% (5%) 


NA 


5% (5%) 


NA 


PDFs/qs /scale 


NA 


7% 


(20%) 


NA 


7% (20%) 


Electroweak K factor 


NA 


2.3% 


(4.5%) 


NA 


2.3% (4.5%) 


Efficiency 


1.0% (2.0%' 


) 1.0% 


(2.0%) 


3.0% (6.0%) 


3.0% (6.0%) 


Scale/Resolution 


1.2% (2.4%; 


) 1.2% 


(2.4%) 


1.2% (12%) 


1.2% (12%) 


Multijet/W+jets background 


NA 


12% 


(26%) 


NA 


< 0.1% 


Total 


5% (6%) 


14% 


(33%) 


6% (14%) 


8% (25%) 



The likelihood of observing a set of n events in N in- 
variant mass bins is given by a product of Poisson prob- 
abilities for each mass bin k: 

C(n\0,P) = l[I^J—. (7 ) 

Li nkl 

According to Bayes' theorem, the posterior probability 
for the parameter 9 given h observed events is 

V(6 | n) = ±£ M (n | 9)P(6), (8) 

where Z is a normalization constant and the marginal- 
ized likelihood £m corresponds to the likelihood after all 
nuisance parameters have been integrated out. This in- 
tegration is performed assuming that the nuisance pa- 
rameters are correlated across all mass bins; Table IIIII 
shows which parameters are taken into account for ei- 
ther or both of the signal and background expectations. 
The prior probability P(9) is chosen to be flat in either 
1/A 2 or 1/A 4 for the CI analysis, and either 1/Mg or 
1/M| for the ADD analysis. These choices are moti- 
vated by the form of Eqs. © and ©. The 95% C.L. 
limit is then obtained by finding the value #i; m satisfy- 
ing §Q llm V{9 | n) d9 = 0.95, where 9 is chosen to be 
1/A 2 , 1/A 4 , 1/M| or 1/Mf. The above calculations 
have been performed with the Bayesian Analysis Toolkit 
(BAT) 61] , which uses a Markov chain Monte Carlo tech- 
nique to integrate over nuisance parameters. 

VII. CONTACT INTERACTIONS ANALYSIS 
AND RESULTS 

To test the consistency between the data and the SM in 
the CI search region (m« > 400 GeV), a likelihood ratio 
test is performed by producing a set of SM-like pseudo- 
experiments and comparing the likelihood ratio between 
the signal+background and pure background hypotheses 
obtained in the data to the results of the pseudoexperi- 
ments. The signal+background likelihood is evaluated at 
the A value that maximizes it. The derived p- value, cor- 
responding to the probability of observing a fluctuation 



in the pseudoexperiments that is at least as signal-like as 
that seen in the data (i.e., with a maximum likelihood ra- 
tio greater than or equal to that obtained in the data) , is 
estimated to be 15% (76%) in the dielectron channel and 
79% (59%) in the dimuon channel for constructive (de- 
structive) interference. These values indicate that there 
is no significant evidence for contact interactions in the 
analyzed data and thus limits are set on the contact in- 
teraction scale A. 

Limits are obtained with the Bayesian method de- 
scribed above. Electroweak corrections are applied to 
both DY and DY+CI samples for consistency, although 
part of the electroweak corrections cannot be computed 
reliably due to the unknown new phenomena represented 
by the contact interaction. This particular choice results 
in slightly more conservative limits. 

The expected 95% C.L. lower limit values on the en- 
ergy scale A are found to be 13.8 ± 1.7 TeV for construc- 
tive interference (A - ) and 10.4 ±1.0 TeV for destructive 
interference (A + ) in the dielectron channel. The cor- 
responding expected limits in the dimuon channel are 
12.7 ± 1.5 TeV and 9.9 ± 1.1 TeV. The quoted uncertain- 
ties correspond to the 68% range of limits surrounding 
the median value (taken to be the expected limit) of all 
limits obtained with a set of pseudoexperiments. Lim- 
its are expected to be stronger in the dielectron channel 
than in the dimuon channel due to the significantly larger 
acceptance for the dielectron selection. 

The observed limits (at 95% C.L.) are A" > 12.1 TeV 
and A + > 9.5 TeV in the dielectron channel for construc- 
tive and destructive interference, respectively. The corre- 
sponding limits in the dimuon channel are A~ > 12.9 TeV 
and A + > 9.6 TeV. These limits are summarized in Ta- 
ble US 

If instead of choosing the prior to be flat in 1/A 2 , it is 
selected to be flat in 1 /A 4 to match the form of the pure 
CI term in Eq. ©, the observed limit in the dielectron 
channel becomes weaker by 0.7 TeV for constructive in- 
terference and 0.4 TeV for destructive interference. The 
corresponding respective shifts to lower values are 1.2 and 
0.6 TeV in the dimuon channel, see Table ITVl 

Finally, a limit is set for the combination of the dielec- 
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TABLE IV. Expected and observed 95% C.L. lower limits 
on the contact interaction energy scale A for the dielectron 
and dimuon channels, as well as for the combination of those 
channels. Results are provided for constructive and destruc- 
tive interference as well as different choices of flat priors: 1/A 2 
and 1/A 4 . 



Channel 


Prior Expected limit [TeV] Observed limit [TeV] 
Constr. Destr. Constr. Destr. 


ee 


1/A 2 
1/A 4 


13.8 
12.5 


10.4 
9.8 


12.1 
11.4 


9.5 
9.1 


/'/' 


1/A 2 
1/A 4 


12.7 
11.6 


9.9 
9.1 


12.9 
11.7 


9.6 
9.0 


ee ± mi 


1/A 2 
1/A 4 


15.0 
13.8 


11.3 
10.5 


13.9 
12.9 


10.2 
9.8 



tron and dimuon channels, assuming lepton universality, 
by computing a combined posterior probability for the 
two channels. The following sources of systematic un- 
certainty are treated as fully correlated between the two 
channels: PDF and as, QCD and electroweak K fac- 
tors, and Z/j* cross section for normalization. All other 
sources are treated as uncorrelated. The resulting com- 
bined limits are A" > 13.9 TeV and A+ > 10.2 TeV for 
the 1/A 2 prior. Table |PVl summarizes all limits for the 
two priors considered in this analysis. 

VIII. LARGE EXTRA DIMENSIONS ANALYSIS 
AND RESULTS 

The search for large extra dimensions is carried out 
similarly to that for contact interactions. A difference 
from the CI analysis is that the DY component present in 
the Sherpa DY+ADD simulated samples is subtracted 
out to compute the net ADD contribution to the total 
event yield. The DY background is modeled with the 
same Pythia DY sample as is used for the CI analy- 
sis. Another difference is that the search is performed 
in only one mass bin with the minimum mass chosen at 
the value giving the strongest expected limit. This opti- 
mization results in a signal region with a minimum mass 
requirement of 1300 GeV as determined from a set of 
pseudoexperiments in each of the dielectron and dimuon 
channels. Table [V] presents the expected and observed 
event yields in the signal region, including the expecta- 
tion for several Ms values in the GRW formalism. 

The consistency between the number of observed 
events in the data and the predicted SM contribution is 
assessed using a set of SM-like pseudoexperiments. Using 
the same likelihood ratio approach as for the CI analysis, 
p- values of 6% and 68% are obtained in the dielectron and 
dimuon channels, respectively. These values indicate that 
there is no significant evidence for large extra dimensions 
and thus limits are set on the scale Ms. The observed 
limits are Ms > 2.73 (2.62) TeV in the dielectron chan- 



TABLE V. Expected and observed number of events with 
m« > 1300 GeV in the dielectron and dimuon channels. 
Yields given for different Ms values correspond to the sum 
of signal and background events, with the signal obtained in 
the GRW formalism. All yields are normalized to the Z peak 
control region. The errors quoted originate from systematic 
uncertainties and limited MC statistics. 



Process 


ee 


/'/' 


DY 


0.89 ±0.21 


0.54 ±0.16 


ti 


< 0.01 


< 0.01 


Diboson 


0.075 ± 0.005 


0.059 ± 0.010 


Multijet/W+jets 


0.16 ±0.20 




Total background 


1.13 ±0.29 


0.60 ±0.16 


M s = 1.5 TeV 


72 ±5 


47 ±9 


Ms = 2.0 TeV 


40.2 ± 2.6 


22 ±4 


M s = 2.5 TeV 


11.7 ±0.9 


6.3 ± 1.1 


M s = 3.0 TeV 


4.2 ±0.4 


2.3 ±0.4 


Data 


2 






TABLE VI. Expected and observed 95% C.L. lower limits on 
Ms in the dielectron and dimuon channels, as well as for the 
combination of those channels without and with the diphoton 
channel in the GRW formalism. Separate results are provided 
for the different choices of flat priors: 1/M| and 1/Mf. 



Channel 


Prior 


Exp. limit [TeV] 


Obs. limit [TeV] 


ee 


1/M S 4 


2.88 


2.73 




1/M S 8 


2.72 


2.62 




1/M| 


2.83 


2.83 




I/Ms 3 


2.61 


2.61 


ee ± [ifi 


1/M| 


3.16 


3.00 




1/M S 8 


2.96 


2.85 


ee + jj.fi -+ 


-77 1/M S 4 


3.43 


3.22 




1/M S 8 


3.27 


3.12 



nel and Ms > 2.83 (2.61) TeV in the dimuon channel at 
95% C.L. with a prior flat in 1/M| (1/Mf). Table EH 
shows these observed limits along with the expected lim- 
its. Limits in the dielectron channel are slightly worse 
than expected due to the larger number of events ob- 
served in the data compared with the SM expectation. 
The dielectron and dimuon channels are combined tak- 
ing correlated systematic uncertainties into account in a 
way identical to the CI analysis. 

A search for large extra dimensions has also been car- 
ried out in the diphoton final state using the data sam- 
ple collected by ATLAS in 2011 0. The results of that 
search are combined with the dilepton results presented 
here with the use of BAT. Correlated sources of system- 
atic uncertainty are treated as follows. The PDF uncer- 
tainty in the SM diphoton and DY production is consid- 
ered to be fully correlated between the ee, and 77 



10 



TABLE VII. Observed 95% C.L. lower limits on M s (in units 
of TeV), including systematic uncertainties, for ADD signal 
in the GRW, Hewett and HLZ formalisms with no K factor 
applied to the signal. Separate results are provided for the 
different choices of flat priors: 1/M| and 1/Mf. 



Channel 


Prior 


GRW Hewett 




















72=3 


n=4 


71=5 


n=6 


71=7 


ee 


1 /Ml 


2.73 


2.44 


3.25 


2.73 


2.47 


2.30 


2.17 




1/Ml 


2.62 


2.48 


2.86 


2.62 


2.49 


2.40 


2.34 


lift 


1/Ml 


2.83 


2.52 


3.36 


2.83 


2.55 


2.38 


2.25 




1/Ml 


2.61 


2.47 


2.85 


2.61 


2.48 


2.40 


2.33 


ee + /i/i 


1/Ml 


3.00 


2.68 


3.57 


3.00 


2.71 


2.52 


2.39 




1/Ml 


2.85 


2.70 


3.11 


2.85 


2.71 


2.62 


2.54 


ee + /j/j 


1/Ml 


3.22 


2.88 


3.83 


3.22 


2.91 


2.71 


2.56 


+ 77 


1/Ml 


3.12 


2.95 


3.40 


3.12 


2.96 


2.86 


2.78 



TABLE VIII. Observed 95% C.L. lower limits on M s (in units 
of TeV), including systematic uncertainties, for ADD signal 
in the GRW, Hewett and HLZ formalisms with K factors of 
1.6 and 1.7 applied to the signal for the dilepton and diphoton 
channels, respectively. Separate results are provided for the 
different choices of flat priors: 1/Ml and 1/Mf. 



Channel 


Prior 


GRW Hewett 






HLZ 














71=3 


n=4 


71=5 


n=6 


n=7 


ee 


1/Ml 


2.95 


2.63 


3.51 


2.95 


2.66 


2.48 


2.34 




1/Ml 


2.82 


2.67 


3.08 


2.82 


2.68 


2.59 


2.52 




1/Ml 


3.07 


2.74 


3.65 


3.07 


2.77 


2.58 


2.44 




1/Ml 


2.82 


2.67 


3.08 


2.82 


2.68 


2.59 


2.52 


ee + /i/i 


1/Ml 


3.27 


2.92 


3.88 


3.27 


2.95 


2.75 


2.60 




1/Ml 


3.09 


2.92 


3.37 


3.09 


2.94 


2.84 


2.76 


ee + /j/j 


1/Ml 


3.51 


3.14 


4.18 


3.51 


3.17 


2.95 


2.79 


+ 77 


1/Ml 


3.39 


3.20 


3.69 


3.39 


3.22 


3.11 


3.02 



channels, whereas the multijet background uncertainty 
is fully correlated between the ee and 77 final states. 
It should be noted that the ee and 77 samples are sta- 
tistically uncorrelated since ee candidates have been ex- 
plicitly removed from the diphoton analysis at the event 
selection stage. The observed and expected combined 
limits are given in Table IVI1 with the most stringent ob- 
served limit obtained for the dilepton-diphoton combina- 
tion: M s > 3.22 (3.12) TeV with a prior flat in 1/M| 
(1/Mf) in the GRW formalism. 

The limits obtained using the GRW formalism have 
been translated into the Hewett and HLZ formalisms us- 
ing Eq. Q with results shown in Table IVII1 Limits are 
also obtained with a K factor applied to the ADD signal 
yield to account for next-to-leading-order QCD correc- 
tions. A constant K factor of 1.6 is applied in the dilep- 
ton channel [62[ and 1.7 in the diphoton channel [63j . 
The dilepton-diphoton combination increases limits by 
approximately 0.2 (0.3) TeV with a prior flat in 1/M| 
(1/Ml), taking QCD corrections into account as shown 
in Table EnH 



IX. CONCLUSIONS 

A search for contact interactions and large extra di- 
mensions has been performed in dielectron and dimuon 
events produced in LHC proton-proton collisions at y/s — 
7 TeV. The data sample corresponds to an integrated lu- 
minosity of 4.9 (5.0) fb _1 of pp collisions in the dielec- 
tron (dimuon) channel recorded with the ATLAS detec- 
tor. No significant deviation from the Standard Model 
is observed in the dilepton mass distributions. Using a 
Bayesian approach with a prior flat in 1/A 2 , as was done 
in most previous searches at hadron colliders, the follow- 
ing 95% C.L. limits are set on the energy scale of contact 
interactions: A" > 12.1 TeV (A+ > 9.5 TeV) in the di- 



electron channel and A' > 12.9 TeV (A+ > 9.6 TeV) 
in the dimuon channel for constructive (destructive) in- 
terference in the left-left isoscalar compositeness model. 
Somewhat weaker limits are obtained with a prior flat 
111 1/A 4 . These limits improve existing bounds on eeqq 
and \x\iqq contact interactions from a single experiment. 
Limits are also set on the scale Ms in the ADD large ex- 
tra dimensions model. Those range from 2.4 to 3.9 TeV 
depending on the choice of model, channel, and prior. 
After combining the dilepton and diphoton searches, the 
limits are in the range from 2.6 to 4.2 TeV. 
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